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ABSTRACT 

We have investigated different scenarios for the origin of the binary millisecond pulsar PSR J1911-5958A 
in NGC6752, the most distant pulsar discovered from the core of a globular cluster to date. The hypothesis 
that it results from a truly primordial binary born in the halo calls for accretion-induced collapse and negligible 
recoil speed at the moment of neutron star formation. Scattering or exchange interactions off cluster stars are 
not consistent with both the observed orbital period and its offset position. We show that a binary system of two 
black holes with (unequal) masses in the range of 3 - 100M Q can live in NGC6752 until present time and can have 
propelled PSR J191 1-5958A into an eccentric peripheral orbit during the last ~ 1 Gyr. 

Subject headings: black hole physics — globular clusters: general — globular cluster: individual (NGC 6752) — 
pulsars: individual (PSR J191 1-5958A) — stars: neutron 
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1. INTRODUCTION 

NGC6752 is a core-collapsed nearby globular cluster (GC) 
containing five millisecond pulsars (MSPs), whose precise pro- 
jected positions in the sky have been recently obtained by 
D' Amico et al. (2002). Three of them reside in the GC core. In 
contrast, the other two pulsars are located in the cluster out- 
skirts. At present, they are the two outermost objects from 
the center of their host cluster known in the catalogue of 41 
pulsars in GCs with accurately determined positions (the third 
ranked object being PSR B2127+11C in M15). The first, PSR 
J1911-5958A (PSR-A according to D' Amico et al. 2002) is a 
binary MSP found at more than 3.3 half-mass radii away from 
the cluster center (corresponding to tpsr-a ^ 8 pc). Its compan- 
ion is a dwarf (probably degenerate) star of mass m* ~ O.2M , 
moving with period P or b = 0.86 days on a circular orbit around 
PSR-A (with eccentricity e < 10" 5 , D'Amico at al. 2001). The 
orbital separation is opsr-a ~ 0.0223 AU, and the binary has a 
gravitational binding energy £j>sr-a ~ 10 47 erg. PSR-A is spin- 
ning with a period of 3.27 ms, has an inferred magnetic field 
of ~ 10 8 Gauss and a characteristic age tpsr-a ~ 15 Gyr; in 
short, it looks like a "canonical" recycled MSP (Lorimer 2001). 
The second, PSR J 1 9 1 1 -6000C, is located at 1 .4 half mass radii 
away from the core and is a MSP similar to PSR-A. However it 
lacks of a companion star and this hints for a scenario in which 
PSR-C has been propelled at its current position by a ionization 
event with a rare high speed star. In this paper we will focus 
our attention on PSR-A as we can constrain its evolution from 
its binary nature. 

PSR-A is 4 times more offset than the neutron star binary 
pulsar PSR B2 127+1 1C in M15 (Anderson et al. 1990), and 
has a much lighter companion. For PSR B2 127+1 1C, Phinney 
& Sigurdsson (1991) explored an ejection scenario in which 
the pulsar exchanged its light donor star with a neutron star 
(NS) flying by. Can we explain the nature of PSR-A in similar 
terms? 

In this paper we consider three possible explanations: (1) 
PSR-A, remnant of a massive star, was the component of a truly 
primordial binary born in the GC halo or expelled from the GC 



core due to the natal kick imparted by the supernova explosion; 
(2) it was ejected from the core via scattering or exchange in- 
teraction^) with cluster stars; or (3) it was propelled into the 
cluster halo by a scattering event involving a black hole binary 
(hereafter [BH+BH]). 

In this letter we will investigate all these three alternatives 
with special consideration to the third, less intuitive, most in- 
triguing hypothesis. This is in part motivated by the recent dis- 
covery that NGC6752 shows an unusual high mass-to-light ra- 
tio >10 in its inner region, corresponding to a mass M core > 
10 4 Mq in the form of low luminosity stellar remnants in its 
core (D'Amico et al. 2002). In addition, X-ray observations 
(Pooley et al. 2001) combined with upper limits on the gas 
content of the GC (D'Amico et al. 2002) seem not to exclude 
that a significant fraction of M core is due to a central (perhaps 
binary) BH of a few hundreds solar masses. 



2. CASE 1: A PRIMORDIAL BINARY ? 

This is the simplest hypothesis among the three. It requires 
orbital stability against dynamical friction, and a low recoil 
speed at the moment of NS formation. Orbital stability depends 
on binary mass, initial apocentric distance r apo and orbital ec- 
centricity, e ec . We computed the dynamical friction timescale 
Tdf for a lie reduction of r apo in the case of a point mass M 
moving in the cluster potential (approximated as a King model) 
under the action of the frictional force (Binney & Tremaine 
1987), considering r apo > tpsr-a- For a centrophobic halo orbit 
(e ec £ 0.5) we find that r df is > 10 n (l .6M /M) yr in NGC6752, 
implying stability of PSR-A and of its progenitor halo binary. 
If the binary has been propelled into a highly eccentric orbit 
(e ec ^0.9) at the time of NS formation, it would remain at the 
periphery for a time t& >7 x 10 8 (1.6M Q /M) yr before return- 
ing to the core. Whenever Tdf < tpsr-a, it is the dynamical 
friction time that sets the conditions on detectability. 

How orbital stability combines with the problem of reten- 
tion of PSR-A in NGC6752? It is known from current proper 
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motion measurements of 13 MSPs belonging to the Galaxy 
(Toscano et al. 1999) that MSPs have relatively large mean 
transverse velocities of 85±13kms _1 and a 3-D peculiar mean 
speed of 130 ± 30kms _1 (Lyne et al. 1998; see also Cordes & 
Chernoff 1997). Theoretical studies aimed at estimating the 
minimum recoil speeds V rec (for those binary MSPs that in- 
voke a common envelope phase preceding the formation of the 
NS; Tauris & Bailes 1996) predict for PSR-A a recoil veloc- 
ity V rec 60km s" 1 , on the basis of its observed orbital period. 
If the cluster escape velocity at the time of the supernova ex- 
plosion was similar to the current value V esc ~ 35 kms" 1 for 
the center of NGC6752, the retention of PSR-A in the cluster 
is indeed problematic. If the supernova took place in the cen- 
tral regions, the binary of PSR-A (of mass >2.5M Q after the 
NS has formed) would have been driven back to the core by 
dynamical friction in < 1 Gyr due to the high orbital eccen- 
tricity. As the supernova blew off in the first Gyr since GC 
formation, the binary would not have survived in the halo until 
present time. PSR-A may, alternatively, descend from a binary 
born in the cluster outskirts with the faintest kick ever observed 
(considering that V esc ~ lOkms" 1 in the halo). This scenario 
calls for the NS of PSR-A being probably formed from the 
accretion-induced collapse of a massive white dwarf (Grindlay 
et al. 1988): in fact in this case collapse is not accompanied by 
severe mass losses and this could minimize the center-of-mass 
runaway impulse speed. In this case, the time is consistent 
with the life-time of PSR-A and its evolution. 

Given the estimated number Af ga i ~ 40,000 of primordial 
MSPs in the Galactic field (e.g. Lorimer 2001) and mass ra- 
tio n < 5 x 10" 4 between all known (~ 200) GCs and the Milky 
Way, it turns out that the overall number of MSPs born in the 
halo from a primordial binary is nearly ~ N ga i(j,f r ~ 3(/ r /0.15), 
where f r corrects for the mass fraction outside radius r in a typ- 
ical GC normalized to its value at r equal to 3 half-mass radii. 
This is a very low number considering that refers to all GCs 
known. 



3. CASE 2: SCATTERING OR EXCHANGE OFF A CLUSTER'S STAR ? 

If PSR-A has been ejected out of the core by a cluster star the 
key issue is whether three-body exchange is favored against in- 
elastic scattering. To address this question, we need to consider 
the binary nature of PSR-A and to distinguish cases in which a 
dynamical interaction occurs before recycling (when the com- 
panion to PSR-A has mass of 1 < m c < 2M Q ) or following it 
(when its mass is m* <~ 0.2M Q ). 

We recall that just at the onset of recycling a "bifurcation" 
orbital period Pbif separates the formation of a converging bi- 
nary (which evolves with decreasing orbital period P or b until 
the mass-losing star becomes degenerate and a short-period bi- 
nary pulsar is formed) from a diverging one (which evolves with 
increasing P or b forming a long period binary pulsar; Tauris & 
Savonije 1999). Given this scenario, PSR-A should result from 
a converging binary. Thus, at the onset of recycling we use an 
upper limit on P or b,i n < Pm ~ 2 days and, as customarily, a mass 
of the donor star ~ 1 - 2 M Q . 

Scattering can occur either before or after recycling as it pre- 
serves, to a high degree, the circularity seen in our binary pul- 
sar. In both cases, a recoil kinetic energy of ~ 3 x 10 46 erg is 
requested for ejecting [NSa+co] (comprising either PSR-A or 
the not yet recycled NS and its companion) into its halo orbit. 



As this must be extracted from the [NSa+co] binding energy, 
recoil at V e j <J V esc imposes an upper limit on the orbital period 
^orb.up which is <; 0.4 days. We infer this value using a collec- 
tion of scattering events extracted from Sigurdsson & Phinney 
(1993): Table 1 gives the binary recoil speed V e j (third column) 
for Porb.up (first line of each case), for a period corresponding 
to the current binary separation (second line), and for Pbif (third 
line). Since PSR-A is a converging system, its period can not 
have ever been shorter than its current value of 0.86 days and 
this excludes the scattering hypothesis. 

Exchange alters the nature of the interacting system and 
makes the new binary rather eccentric and with a longer orbital 
period. Thus due to the prohibitively long circularization time 
between compact stars, exchange has to occur before recycling. 
We consider the case (see Table 1) of a main sequence binary 
[MS+ms] (MS heavier than ms) interacting with the NSa fly- 
ing by (given the similar masses of MS and NS A , the results 
are almost unaffected when swapping their role). We find that 
the target binary needs to be very tightly bound to produce the 
new system [NSa+co]; after the exchange [NSa+co] has an or- 
bital period suspiciously close to the currently observed value 
considering that recycling has not occurred yet; characteristic 
times ~ 1.8 x 10 10 years for exchange (see Table 1; first line of 
both exchange cases) are also found uncomfortably long in this 
case. 



4. CASE 3: SCATTERING OFF A [BH+BH] ? 

The dynamics of BHs in GCs involves repeated exchange in- 
teractions and ejections (Sigurdsson & Hernquist 1993; Porte- 
gies Zwart & McMillan 2000). This cosmic dance occurred 
early in the cluster lifetime among the rich population (of few 
hundreds) of BHs, born from the most massive stars: As they 
sink toward the center by dynamical friction (on a time-scale 
<J 10 6 yr), BH-[BH+BH] interactions become overwhelmingly 
important. Binaries harden progressively in virtue of these en- 
counters and eventually leave the cluster due to recoil. As a re- 
sult, many clusters lose their original BHs on a time scale of ~ 
Gyr. However, dynamical studies (Portegies Zwart & McMillan 
2000; Miller & Hamilton 2002a) show that some clusters could 
still retain a binary. We here consider the fascinating possibil- 
ity that NGC6752 belongs to this category and try to depict a 
scenario for the interaction of PSR-A with a [BH+BH] system. 

Two separate questions arise: (i) For how long can the 
[BH+BH] binary that survived ejection avoid coalescence by 
emission of gravitational waves (GWs)? (if) Can the [BH+BH] 
transfer enough kinetic energy to propel PSR-A into the cluster 
halo and what is the probability of this event? 



4. 1 . GW Emission and Recoil in the [BH+BH] 

Here we allow BH masses in binaries to have a relatively 
wide spectrum from 3 M Q to 100M Q . The BHs may result from 
the direct collapse of massive progenitor stars (Fryer 1999). But 
the massive >50M Q -BHs could be the outcome of exchanges 
and mergers of many ordinary (~ 10M Q ) BHs; in particular 
binary-binary encounters can enhance their formation (Miller 
& Hamilton 2002b). As noticed by Miller & Hamilton (2002a), 
unlike the less massive BHs in binaries which can be easily 
flung out by recoil before they can merge, the more massive 
ones could have enough inertia to remain bound to the cluster. 
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Thus, the binary that avoided ejection has likely unequal masses 
and could harbor the heaviest BH that NGC6752 ever had. 

If we indicate with Mbh and w B h the masses of the two BHs 
in the binary, respectively, and with M Bli t and /i BH their to- 
tal and reduced mass, our hypothetical [BH+BH] binary (with 
typical eccentricity of 0.7) will avoid coalescence for a time 
t"gw > 10 Gyr due to emission of GWs if its separation a ex- 
ceeds a critical value ogw, where 



«gw ~ 0.4 



/ Af] 



BH,t 



V 100M, 



1/2 



/ Mbh 



V 1OM 10 10 yr 



1/4 



AU. (1) 



Recoil, instead, imposes a different condition based upon the 
request that the BH separation is sufficiently large to avoid ejec- 
tion: the [BH+BH] binary that remains bound to NGC6752 
should have a > a lec , where 



0.01 (-!^- 

Vi0M Q 



m 100M Q 35kms 
10M o M BH .t V^T 



AU. (2) 



Equation (g) is derived knowing that binaries interacting with 
much lighter intruders (of mass m) acquire a kinetic energy 
(due to recoil) of the order of ~ 0.2(m/MBH.t) 2 £B, where £b = 
GMBH'WBH/2flrec is the binding energy of the [BH+BH] that 
typically produces a recoil at the escape speed V esc ~ 35kms _1 
from the cluster core (see Miller & Hamilton 2002a; Hills 1992; 
Quinlan 1996). In Figure 1 we show these critical separations; 
when flow > flrec coalescence is the correct criterion to use. 

After the last encounter with a background BH, the evolution 
of [BH+BH] continues; what comes into play is its interaction 
with the bath of cluster stars, of mass (m) , line of sight velocity 
(v), and density n*. If £* ~ 2iraGMsu x / (v 2 ) is the three-body 
interaction cross section, then the (still very approximate) num- 
ber of stars impinging onto the target binary over a time t ~ 10 10 
yr is A/"* ~ (v)t and gives 



KT <l U$f n * Mbh ' 1 6kmS ' a 



t 

(v) flew 10 To yr 



(3) 



for NGC6752. Each single star carries off an energy ~ 
0.2((m) /MsH.d^B so that the [BH+BH] hardens in the course 
of time reaching a = agw from an initially much larger sep- 
aration a acc ~ 2GM Blit / (v 2 ) (taken to be comparable to the 
accretion radius) after experiencing a number of scatterings 
A/gw ~ 5[M B H,t/ (»i)]ln(a aC c/flGw) that exceeds A/"* 



GW 



^5-10 



/ M BH .t O-4M \ / q acc 0.4AU \ 
V1OOM (m) J n ^80AU a GW / 



(4) 



This secures a [BH+BH] lifetime comparable to that of 
NGC6752, or slightly shorter. 
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Fig. 1. — [BH+BH] binary separation a against itlbh / Mq for Mbh = 
50Mq ; ogw (heavv solid line), and a scat (heavy dashed line) are com- 
puted from eqs. 0, and |^. a rcc , from eq. g is computed for light 
m = IOMq BH intruders and for m = otbh- The dotted line gives the 
tidal disruption radius for [NSa+co]. The shaded area indicates the 
values of a for ejection of PSR-A during the last Gyr, compatible with 
its detection. 



4.2. PSR-A Ejection 

Let us now come to question («). In the ejection hypothe- 
ses, the binary [NSa+co] must have been propelled into an ec- 
centric peripheral orbit at a speed V e j < V esc . The apocenter 
then reduces due to dynamical friction but its e-folding time Tdf 
>7 x 10 8 (1 .6M Q /M) yr is sufficiently long to permit the detec- 
tion of PSR-A during the pulsar life-time ~ tpsr-a- 

In the scattering off [BH+BH], the light binary [NSa+co] is 
viewed as a single point mass. In this case, [NSa+co] has to 
avoid break up by the BH tidal field and this imposes a con- 
straint on the distance of closest approach, r m ; n > r t id where 



r tid ~ 0.2 



M 



BH.t 



l.6Mr. 



IOOMq m PSR _ A +m c 



1/3 

AU. 



(5) 



In the inelastic scattering off the hard [BH+BH], most of the 
released binding energy Ae^ ~ 0.2[(mpsR-A + fw c )/MBH,t] £b of 
the binary goes into kinetic energy of the propelled [NSa+co]. 
Thus for an ejection velocity Vej close to Vesc, the [BH+BH] 
should have a binary separation smaller than 



^scat 



V 2 



1.3 



V 10M, 



35kms 

Vei 



AU. (6) 



Since scattering leads to an effective energy exchange only 
when the intruder ([NSa+co] in this case) is approaching the 
binary within a distance r comparable to a few orbital separa- 
tion radii a (beyond which energy exchange is exponentially 
small), conditions (|||§) combined with (|l|) or (Q) give a plausi- 
ble range of binary separations for which injection of [NSa+co] 



4 



into the cluster halo is possible and this is illustrated in Fig. 1 
(shaded area). 

If the flux of stars on the [BH+BH] is of thousands interac- 
tions over the cluster lifetime (see eq. || and Q), about 1/30 of 
these encounters should have involved a NS, which gives a rate 
1Z„ S ~ 10 Gyr -1 (taking as face value 3,000 encounters), but it 
can be higher (lower) if the [BH+BH] mass is higher (lower). 
If — 5% of the NSs is in the state of a MSP, then the rate for 
MSP ejection is IZmsp ~ 0.5 Gyr -1 , implying a probability of 
detection of about ~ 50% considering a life-time for the cur- 
rent orbit of ~ 1 Gyr. We note that 1Z„ S is calculated assuming 
a Salpeter-like initial mass function (giving ~ 1% of NSs) and 
a mass segregation enhancement of ~ 6 in the density (Sigurds- 
son & Phinney 1995). If the IMF is flatter (as suggested by the 
observations of D' Amico et al. 2002), this rate could be higher, 
enhancing the likelihood of this scenario. Once ejected into an 
halo orbit, [NSa+co] spends only ~ 1/50 of its life in the dense 
GC core and thus has negligible probability to undergo further 
interactions (see Table 1). 

5. DISCUSSION 

The unusual position of PSR-A in NGC6752 is puzzling. We 
have shown that PSR-A may result from a primordial binary 
born in the halo, but this requires an exceptionally weak kick 
or accretion-induced collapse at the moment of NS formation. 



A scattering or exchange event off a star that propelled PSR- 
A into the cluster halo is inconsistent with the conservation 
of energy and linear momentum. Here, we propose a dynam- 
ical interaction with a more massive system, a [BH+BH] bi- 
nary in a mass range of ~ 3- 100M Q and separation of ~ 0.3-2 
AU. An unequal mass binary with a heavy >100M Q BH and a 
lighter companion ~ 10M Q is preferred, as it can form in a GC 
(Miller & Hamilton 2002a, 2002b) and can provide the right 
statistics for explaining PSR-A. Light binaries are likely to be 
ejected (Portegies Zwart & McMillan 2000), while much heav- 
ier >2OOM ones (with equal mass) are difficult to form and 
their presence can be in conflict with X-ray observations if the 
gas content in NGC6752 is confirmed at the density reported 
in D'Amico et al. (2002); in addition, the dynamical effects of 

such a heavy binary on cluster stars could also become severe. 

l 

Whether a GC retains its "last" BH binary(ies) or not is un- 
predictable, as BH ejection depends on a number of rare ran- 
dom events. It is thus an observational challenge to unveil in 
the heart of a cluster such an exotic dark binary. NGC6752 
seems a primary candidate for this search. 
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MM02C7 1842-001. 
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Table 1 

Dynamical interactions for PSR-A ejection 



forb (I V'ej 

(d) (AU) (kms^ 1 ) 


T 

(yr) 


a ' Porb 

(AU) (d) 


scattering (b) (1.4,1.4MS) 


+ (0.6ms) -> (1.4, 1.4MS) + (0.6ms) 


0.089 0.0055 33 
0.713 0.0220 16 
2.018 0.0440 12 


4.6 x 10 s 
5.8 x 10 8 
2.1 x 10 8 


0.0054 0.086 
0.0215 0.689 
0.0430 1.948 


scattering lb, (1.4, 1MS)h 


r(1.4MS)- 


»(1.4,1MS) + (1.4MS) 


0.361 0.0135 34 
0.752 0.0220 26 
2.019 0.0425 19 


2.8 x 10" 
1.4 x 10 9 
5.1 x 10 8 


0.0131 0.344 
0.0213 0.715 
0.0411 1.920 


scattering (a, (l .4, 0.5ms) + (0.7ms) - 


» (1.4, 0.5ms) + (0.7ms) 


0.067 0.0040 33 
0.868 0.0220 14 
2.010 0.0385 10 


1.6 x 10 1U 
1.2 x 10 9 
5.2 x 10 8 


0.0039 0.065 
0.0216 0.843 
0.0378 1.952 


scattering (a, (l .4, 0. 14ms) H 


-(1.4MS) - 


>(1.4,0.14ms) + (1.4MS) 


0.231 0.0085 33 
0.962 0.0220 21 
2.013 0.0360 16 


2.9 x 10 1U 
6.9 x 10 9 
3.3 x 10 9 


0.0078 0.205 
0.0203 0.853 
0.0332 1.786 


exchange (b, (0.7MS, 0.5ms) + (1.4) - 


>(0.7MS, 1.4) + (0.5ms) 


0.382 0.0118 34 
1.029 0.0229 24 
2.029 0.0360 19 


6.0 x 10 9 
2.2 x 10 9 
1.2 x 10 9 


0.0330 1.424 
0.0640 3.836 
0.1006 7.563 


exchange (b, (1.4MS, 0.6ms) + (1.4) - 


>(1.4MS,1.4) + (0.6ms) 


0.382 0.0130 33 
0.865 0.0224 25 
2.029 0.0396 19 


6.0 x 10 9 
2.7 x 10 9 
1.2 x 10 9 


0.0318 1.2389 
0.0548 2.805 
0.0968 6.580 



NOTE. — Scattering and exchanges from Sigurdsson & Phin- 
ney (1993); heavy (lighter) main sequence stars are denoted with 
MS (ms). When indication is absent the star identify with NSa of 
PSR-A. Process labeled with b (or a) occurs before (or after) re- 
cycling. Vej, a', and P' mb denote the ejection velocity of [NSa+co], 
the binary separation and orbital period, after the encounter. T is 
the interaction time scale in NGC6752 for n = 10 6 pc~ 3 (D' Amico 
et al. 2002); for exchanges with incoming NSs the value of T 
given in the table needs to be multiplied by a factor ~ 30, equal 
to the expected ratio between the stellar to NS density. 



